Interfaces between highly incompatible polymers 

of different stiffness: 
Monte Carlo simulations and self-consistent field 

calculations 

M. and A. Werner^ 

^ Institut fiir Physik, Johannes Gutenberg Universitat 
D-55099 Mainz, Germany 
■^Department of Physics, Box 351560, University of Washington, 
Seattle, Washington 98195-1560 

February 1, 2008 

Abstract 

We investigate interfacial properties between two highly incompatible 
polymers of different stiffness. The extensive Monte Carlo simulations 
of the binary polymer melt yield detailed interfacial profiles and the in- 
terfacial tension via an analysis of capillary fluctuations. We extract an 
effective Flory-Huggins parameter from the simulations, which is used in 
self-consistent field calculations. These take due account of the chain ar- 
chitecture via a partial enumeration of the single chain partition function, 
using chain conformations obtained by Monte Carlo simulations of the 
pure phases. The agreement between the simulations and self-consistent 
field calculations is almost quantitative, however we find deviations from 
the predictions of the Gaussian chain model for high incompatibilities 
or large stiffness. The interfacial width at very high incompatibilities is 
smaller than the prediction of the Gaussian chain model, and decreases 
upon increasing the statistical segment length of the semi-flexible compo- 
nent. 



1 Introduction 

Melt blending of polymers has proven useful in designing new composite ma- 
terials with improved application properties. In many practical situations the 
constituents of the blend are characterized by some degree of structural asymme- 
try. For example, a flexible component might contribute to a higher resistance to 
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fracture, while blending it with a stiffer polymer can increase the tensile strength 
of the material. Since the entropy of mixing in polymeric systems decreases with 
increasing degree of polymerization, a small unfavorable mismatch in enthalpic 
interactions, entropic packing effects or the combination of both, generally leads 
to materials which are not homogeneous on mesoscopic scales, but rather fine 
dispersions of one component in another. Therefore properties of interfaces be- 
tween unmixed phases are crucial in controlling the application properties of 
composites jll and have found abiding experimental interest [0, |[ |, ||. 

Recently, the bulk phase behavior and surface propertiesP| of polyolcfins|0, 
^ with varying microstructure has attracted considerable experimental and the- 
oretical interest. These mixtures are often modeled ^ |l^ as blends of poly- 
mers with different bending rigidities, the less branched polymer corresponding 
to the more flexible component. For pure hard core interactions, field theoretical 
calculations by Fredrickson, Liu and Bates polymer reference interaction site 
model (P-RISM) computations by Singh and SchweizerJTo|, lattice cluster the- 
ories by Freed and Dudowiczjll and Monte Carlo simulations]!^ find a small 
positive contribution to the Flory-Huggins parameter x- Monte Carlo simu- 
lations which include a repulsion between unlike species reveal an additional 
increase of the effective Flory-Huggins parameter with chain stiffness, because 
a back folding of chains becomes less probable with increasing stiffness and the 
number of intermolecular contacts increases |12| r espectively. Qualitatively simi- 
lar effects were found analytically in P-RISM|10| and lattice clusterjll theories. 

In spite of their ubiquitous occurrence, interfacial properties in asymmet- 
ric blends have attracted comparably little interest. When entropic packing 
contributions to the Flory-Huggins parameter x a-i'e small and composition fluc- 
tuations are negligible, the self-consistent field theory is expected to yield an 
adequate quantitative description. Helfand and Sapsep^ extended the self- 
consistent field theory to Gaussian chains with different statistical segment 
lengths. In the limit of infinite long Gaussian chains and strong segregation, 
they obtained analytical expressions for the interfacial width w and the inter- 
facial tension a. Both increase upon increasing the statistical segment length 
of one component, leaving x s-i^d the architecture of the other component unal- 
tered. 

However, there are other models, that incorporate structural disparities on 
the monomer level. Freed and coworkers model monomers as clusters of various 
shape on a lattice IT^ and have explored corrections to the energy of mixing and 
entropic contributions to the Flory-Huggins parameter. 

Stiffness disparities have also been investigated using the worm-like chain 
model [|l5| , which captures the crossover between rod- like behavior on small 
length scales and Gaussian statistics on length scales much larger than the 
persistence length. Morse and Fredrickson[^ extended the self-consistent field 
calculation to a symmetric blend of worm-like chains. For vanishing bending 
rigidity k they reproduced the Gaussian chain result. In the limit of high bend- 
ing rigidities and strong segregation {nx ^ 1), however, they found that the 
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width w of the monomer density profile can be considerably smaller than for a 
Gaussian chain with the same long distance behavior. At large kx, increasing 
the statistical segment length even leads to a decrease of the interfacial width in 
qualitative contrast to the Gaussian chain result. They also observed that the 
width of the bond orientation profile is of the order of the persistence length, 
which is much larger than w in that limit. Thus the interfacial width w and 
the persistence length constitute two independent length scales of the interfa- 
cial profiles. A reduction of the interfacial width in the case of small bending 
rigidities was obtained numerically by Schmid and Miillcrjl^. They noted that 
the local structure might become important if its length scale is comparable to 
the interfacial width; a situation which occurs at rather large incompatibilities. 

In the present study we extend our Monte Carlo studies of structural 
asymmetric blends to the investigation of interfacial properties between well 
segregated phases of flexible and semi-flexible polymers. We consider rather 
small bending rigidities of the semi- flexible component, so that the long distance 
behavior of both species is Gaussian. However, we chose the incompatibility x 
high enough, such that the interfacial width and the persistence length are 
comparable for the higher bending rigidities. The Monte Carlo simulations 
highlight the architectural influences and give a detailed picture of interfaces 
between structural asymmetric polymers. They yield density and orientation 
profiles for bonds and chains as a whole. Extracting an effective Flory-Huggins 
parameter x from the simulation data, we compare our Monte Carlo results to 
self-consistent field calculations which take due account of the chain architecture 
via a partial enumeration procedure p8|, p^ , pc| ], and to Gaussian chain results. 
Therefore we can assess the importance of the level of coarse graining on the 
interfacial properties. 

Our paper is organized as follows: In the next section we describe our poly- 
mer model, especially the dependence of single chain properties on the stiffness. 
We comment on some computational aspects of the Monte Carlo simulations 
and describe the measurement of the interfacial tension. We also introduce the 
salient features of our self-consistent field calculations for arbitrary molecular 
architecture. In the following, we present our simulational results and compare 
them to the self-consistent field calculations. We close with a brief discussion 
of our findings and an outlook on future work. 

2 Model and technical details 

2.1 Bond fluctuation model and single chain properties 

In the framework of our coarse grained lattice model, a small number of chem- 
ical repeat units, say 3-5, is mapped onto a lattice monomer, such that the 
relevant features - chain connectivity and excluded volume interaction between 
monomeric units - are retained. We use the three dimensional bond fluctua- 
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tion model (BFM)||2l|, which has fomid widely spread apphcation in computer 
simulations, because it combines the computational efficiency of lattice models 
with a rather faithful approximation of continuous space properties. Each effec- 
tive monomer blocks a cube of 8 neighboring sites from further occupancy on 
a simple cubic lattice. Due to the extended monomer size, the model captures 
some nontrivial packing effects. We consider a blend of 71^ flexible polymers 
of length and semi-flexible B-polymers comprising Ng monomers in a 
volume V. At a total monomer density $0 = (A^a^a + ^^8^3)/^ = 0.5/8, the 
model reproduces many properties of a dense polymeric melt. We use chain 
lengths N — Na = Nb — 32 and 64, which correspond to a degree of polymer- 
ization of the order 120 and 240 in more chemically realistic polymer models. 
Monomers are connected via one of 108 bond vectors with lengths 2, 3 
or -s/lO, where here, and henceforth, all lengths are measured in units of the 
lattice spacing. The large number of bond vectors permits 87 different bond 
angles. 

The persistence length of the semi-flexible B-polymers is tuned by imposing 
an intermolecular potential, which favors straight bond angles. We use a partic- 
ular simple choice [T2|: E{0) = fksT cos{6) where 6 denotes the complementary 
angle to two successive bonds. Previous Monte Carlo simulations |p2[ of the 
bulk thermodynamics for = 32 and / = 1.0 revealed a purely entropic Flory- 
Huggins parameter Ax — 0.0018(2) for the athermal blend. This small value is 
in good quantitative agreement with theories j9[ |l^. These packing effects result 
in a slight increase of the osmotic pressure with the bending energy, which gives 
rise to a monomer density difference of about 1% between the coexisting phases. 

Since Ax <C 2/A^ — 0.0625 for this combination of chain length and stiffness 
disparity, we introduce an additional enthalpic repulsion to induce phase sepa- 
ration. For simplicity, these thermal interactions are modeled as a square well 
potential comprising all 54 neighbor sites up to a distance -\/6- The contact of 
monomers of the same species lowers the energy by e/csT, whereas the contact of 
different monomers increases the energy by the same amount. A finite size scal- 
ing analysis yields accurate estimates for the critical point of the binary blends 
{N = 32): ec = 0.01442(6), = 0.5 and = 0.0127(1), 0ac = 0.516(10) for 
/ = and ijl^, respectively. In the present study we chose e = 0.05 which 
corresponds roughly to x ~ 0.27. This value is much higher than typical values 
for polyolefin blends Q. Our results correspond to rather strongly immiscible 
blends (e.g. interfaces between polystyrene (PS) and polyvinylpropylene PVP 

i). 

The conformational data for A^ = 32 and e — 0.05 as a function of the 
bending energy / are presented in Fig. |] and Table |l|. The inset shows the 
growth of the chain extension upon increasing /. The ratio between the square 
end-to-end distance and the square radius of gyration it!^ remains very close 
to the Gaussian value 6 (within 5% even for f ~ 2). Also the small wave vector 
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regime of the single chain structure function 
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is well describable by a Debye function SD{q)/N = 2[exp{-q^Rl)-l+q^Rl]/{q^Rl)^^ 
for q < 0.3. Thus the long range behavior of our chains is characterized by Gaus- 
sian statistics for all values of the bending energy / studied and we define the 
statistical segment length b according — /{N ~ 1). Note that the statistical 
segment length grows from 3.06 for / = to 4.63 for f — 2. This asymmetry in 
the statistical segment length is of similar magnitude as in polyolefin blends |Q. 

However, for length scales of the order of the statistical segment length, we 
find deviations from the Gaussian behavior. The plateau q^S{q) — 12/6^ for 
large q in the Kratky Porod plot is only observed for flexible chains (/ = 0) 
and yields a slightly higher estimate for the statistical segment length b = 
3.4. For the semi-flexible chains the slope of q'^S{q) in the range 0.3 < g < 1 
increases upon increasing the bending energy. Defining an effective bending 
rigidity of an equivalent worm-like chain k = R^ /2{b'^){N — 1) ((6^): mean 
squared bond length), k grows from 0.68 to 1.57 [kx = 0.18- ••0.47) upon 
increasing the bending energy /. For wave vectors qRg ~ 27ri?g/2w (denoted 
by the arrows in the Fig. |l|), where w = 3.4 corresponds roughly to the width 
of the monomer density profile in the self-consistent field (SCF) calculations, 
we find deviations from the Gaussian behavior for higher bending energies and 
anticipate corrections to the predictions of the Gaussian model. 

2.2 Local fluid structure and efl'ective Flory-Huggins pa- 
rameter X 

In order to compare our simulational results to self-consistent field (SCF) calcu- 
lations, which cannot account for the local fluid structure of our model, we have 
to identify an effective Flory-Huggins parameter x- For the bulk behavior in 
the one phase region this has been discussed in ref . ^ : We define a dimen- 
sionless monomer density 4)a{b) as the ratio between the local number density 
of A(B)-monomers and the total monomer density ^q. Then, the density of 
intermolecular contacts uab takes the form: 



where gu denotes the IJ interchain correlation function, which is normalized 
such that gij{r ~f oo) = 1. The integration is extended over the spatial exten- 
sion of the square well potential and z/j corresponds to the effective coordina- 
tion number of the Flory-Huggins lattice. If the coupling between chain confor- 
mations and effective monomer repulsion is negligible, only the mterarolecular 
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energy drives the phase separation. In this case (as we shall see in the next 
subsection), the x parameter takes the form: x — ^{zaa + 2zas + zbb) /'^.kBT , 
where zij denote the coordination numbers obtained from the znteTinolccular 
pair-correlation functions. At the critical temperatures the coordination num- 
bers have been measured in the simulations at composition (j)A = 4>b = l/2p^. 
From its very definition the Flory-Huggins parameter x accurately describes 
the intermolecular interaction energy, and it agrees nicely with values obtained 
from the semi-grandcanonical equation of state and the estimate from the long 
wavelength behavior of the collective structure factor. It also yields estimates of 
the critical temperature, which agree with the Monte Carlo results up to 1/VN 
corrections due to composition fluctuations[p3|. 

In the pure system, the intermolecular coordination number of the flexible 
component is lower than the corresponding value for the semi- flexible chains fl^ . 
The number of intramolecular contacts |Q is higher for the flexible chains. 
Therefore, the ^-parameter grows upon increasing /[|T2[. Due to the larger chain 
extension for the semi- flexible component, the correlation hole has a larger spa- 
tial extent, but is more shallow. The intermolecular pair correlation function is 
presented in the inset of Fig. ^ Due to the extended monomer size g{r) vanishes 
for distances r < 2. At short distances, the presence of single site vacancies in- 
troduces local packing effects, which gives rise to several neighbor shells in the 
fluid. The extended structure of the polymer manifests itself in a reduction of 
contacts with other chains on the length scale of the end-to-end distance. On 
short distances, the intermolecular pair-correlation function for the stiffer chains 
is larger than for the flexible ones. For flexible chains it is possible to separate 
the monomeric packing effect from the polymeric correlation hole by dividing 
g{r) by its monomeric equivalent |23[| , which exhibits only packing effects. The 
ratio gir)/ gN=i{f) presents the conditional probability of flnding a monomer of 
a different chain at a distance r, if there would be one in the monomer system. 
This ratio, presented in Fig. |^, is a rather smooth function, indicating, that the 
chain connectivity hardly affects the monomeric packing. If the correlation hole 
would be characterized by a single length scale, i.e. the end-to-end distance R 
in the Gaussian chain model, one expects a scaling behavior of the form: 



Such a scaling plot is shown in Fig. The data collapse well for the different 
bending rigidities at large distances, whereas there are deviations for small dis- 
tances. This is a further indication, that the chain structure is characterized 
by two independent length scales, the end-to-end distance and the persistence 
length. 

In the well segregated regime (far below the critical temperature), it is 
very difficult to measure the AB intermolecular correlation function in the 
bulk. Therefore, unlike ref. [l^ , we make an additional ad-hoc assumption: 
zab — {zaa + zbb)/'2- For symmetric blends near the critical point P-RISM 
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calculations [p5[ predict that deviations from this behavior die out with grow- 
ing chain length like \/^/N. However, the validity of this random-packing like 
assumption for highly incompatible structural asymmetric blends is not obvious. 

We explore the interfacial structure by simulating a system in a L x L x 
2L geometry with L = 64 and periodic boundary conditions in the canonical 
ensemble. The system contains two interfaces parallel to the xy plane. The chain 
conformations are generated via local monomer displacements and slithering 
snake moves, which are applied at a ratio 1:3 (except for / = 0, where only 
local monomer displacements were employed). The systems were equilibrated 
over 125,000 attempted local moves per monomer (AMM) and 375,000 slithering 
snake tries per chain (SS). Every 12,500 AMM and 37,500 SS movements a 
configuration was stored for detailed analysis, at least 898 configuration were 
generated. We use a trivial parallelization strategy on a CRAY T3E, running 
typically 8 or 32 configurations in parallel. 

Profiles across the interface are measured according to the following proce- 
dures: "Apparent" profiles are obtained by locating the instantaneous position 
of the interface across the whole lateral system extension in each snapshot and 
averaging over profiles with respect to the instantaneous, but laterally averaged 
midpoint. These profiles exhibit a system size dependent broadening due to cap- 
illary fluctuations, which is not accounted for in the SCF calculations. To avoid 
this broadening, we define "reduced" profiles by laterally dividing the system 
into subsystems of size B x B. We choose B — 16. One could reduce the effect 
of capillary fluctuations further by chosing a smaller block size B, however, one 
should take care not to cut off "intrinsic" fluctuations j2^ . Since on the scale 
B fluctuations are still reasonably described by a Helfrich Hamiltonian|27[ (see 
below), our block size B is larger than the length scale of "intrinsic" fluctu- 
ations. This is consistent with Semenov's|2^ estimate for the corresponding 
length scale Lcutoff = t^w Ri 10 < B. Thus this averaging procedure reduces the 
influence of capillary fluctuations, but does not eliminate it completely [p8|. 

The presence of an interface gives rise to a spatial dependence of the local 
monomer densities and chain conformations, which in turn is reflected in the 
intermolecular pair correlation functions. In Fig. ^ we present the reduced 
proflles of the intermolecular and intramolecular coordination numbers as a 
function of the distance from the center of the interface for the bending energies 
f — and 2. The individual coordination numbers exhibit a considerable spatial 
dependence; this is however partially due to the spatial range of interactions 
and the remaining capillary fluctuations |2^ . To illustrate the effect we plot the 
apparent and reduced profile of the AB intermolecular coordination number. 
The value at the center of the interface increases upon reducing B; the intrinsic 
value can not be estimated from these data with high precision. However, the 
average value of the "reduced" profile is close to zab = {zaa + zbb)/'2., the 
value used in the SCF calculations. 

The total number of intermolecular contacts z™*'^'' — {n'^AA '' + ^bb^ + 
u'ab^) /^o{4'A + 4>bY is much less sensitive to the intrinsic (local) profiles and 
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shows a gradual transition between the corresponding bulk values, with a re- 
duction at the center of the interface Q of about 8%. This spatial dependence 
of the effective Flory Huggins parameter is neglected. Interestingly, the sum 
of all contacts (both intermolecular and intramolecular) is largely indepen- 
dent of the stiffness or the distance from the interface, i.e. the bending energy 
or the unfavorable interactions at the interface causes the chains to rearrange 
(e.g. exchange unfavorable interchain contacts by energetic favorable intrachain 
contacts) but hardly affect the structure of the underlying monomer fluid. 

Therefore, the local fluid structure is dominated by the packing constraints 
and the excluded volume interactions. The chain connectivity, bending energies, 
and the thermal interactions are of minor importance for the monomer fluid. 
The chain conformations are strongly influenced by the bending energies but 
depend only slightly on the thermal interactions. The Flory-Huggins parameter 
is determined by the thermal interactions and also depends on the bending 
energies via the correlation hole effect. The disparity in the packing behavior 
of the flexible and the stiff polymers is of minor importance for x for the chain 
lengths studied. 



2.3 Measuring the surface tension via the capillary fluc- 
tuation spectrum 

Due to the stiffness disparity between the species, straightforward application of 
semi-grandcanonical identity changes between different polymer types are rather 
inefficient (note that the efficiency drops by about 3 orders of magnitude [p^ 
upon increasing / from to 1 for N — 32) and therefore limited to small chain 
length and stiffness. Measurement of the interfacial tension via the reweighting 
of the composition distribution, which has been successfully applied to structural 
symmetric blends (/ — 0), is therefore difficult. In principle, the interfacial ten- 
sion can be determined via the anisotropy of the pressure tensor. This method 
has been successfully applied in off-lattice simulations , but the generalization 
to lattice models is difficult [Q. However, the spectrum of capillary fluctuations 
offers an alternative |3^ for measuring the interfacial tension; a method which 
does not rely on identity switches. Let u{x,y) be the local interfacial position. 
Then the free energy cost for deviations from a flat planar interface is given by 
the Helfrich expression P7|: 



n 



dxdy -^{^uf + 



(4) 



where higher order gradient terms are neglected. In our simulation, we define 
local X- and y-averaged interface positions by minimizing the quantity 
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for the x-averaged position u(i/) and a similar expression for the y-averaged 
one. This averaged interfacial position is Fourier decomposed according to: 
"(y) = ^ + Y.i=o a('7i)cos(g;2/) + b{qi) sm{qiy) with qi = 2itI/L. The Helfrich 
Hamihonian predicts that the Fourier components a{qi) and h{qi) are Gaussian 
distributed with a width 

In Fig. Q we present the distribution of the Fourier components for two dif- 
ferent bending energies / = 0, 2 and the 4 smaUest wave vectors q. This fong 
wavelength part of the fluctuation spectrum is well described by the quadratic 
Helfrich expression. The straight line marks the expected Gaussian distribution 
for the Fourier amplitudes, to which the simulation data comply. The inverse 
width of the distribution determines the interfacial tension. The extracted value 
for the symmetric blend agrees with the independent measurement obtained via 
the reweighting schemejQ. (The latter scheme measures the interfacial free 
energy via the ratio of the probability for finding the system in a homogeneous 
bulk state or a configuration comprising two interfaces. ) To estimate the er- 
rors of measuring the interfacial tension via the capillary fluctuation spectrum, 
it would be desirable to increase the lateral system size. However, the error 
in extrapolating the simulation data to g ^ is smaller than 7%. Thus the 
analysis of the capillary fluctuation spectrum is an efflcient alternative for mea- 
suring interfacial tensions in structurally asymmetric systems; the results are 
compared to the predictions of the SCF calculations in Sec. Ilia. 



2.4 Self-consistent field calculations 

The mean field approach is similar to Helfandp^, p^, Noolandip5||, and ShuUQ, 
except for the treatment of the chain architecture The partition function 
of a binary polymer blend has the general formJSTt: 

(7) 

where the functional integrals I?[r] sum over all polymer conformations and V[r] 
denotes the probability distribution characterizing the noninteracting, single 
chain conformations. £ represents a segmental interaction free energy, and the 
dimensionless monomer density takes the form[^: 



a—l i/ 



where the sum runs over all monomers in the A-polymer a. A similar expression 
holds for (pBir). 
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The segment free energy £ comprises two contributions: a free volume part 
arising from hard core interactions and an energetic term from the thermal in- 
teractions. Since the melt is nearly incompressible, we approximate the free 
volume part by a simple quadratic expression introduced by Helfand|Q, which 
reproduces the relative reduction of the total monomer density by about 4%|jl^. 
However the difference of the bulk densities of the coexisting phases has a differ- 
ent sign in the simulations and than in the SCF calculations. In the simulations 
the higher osmotic pressure of the semi-flexible component results in a slightly 
lower bulk density of the semi-flexible component in the simulations, an effect 
neglected in the SCF calculations. Moreover, in the SCF calculations the more 
negative intermolecular energy density (see below) of the B component results in 
a slightly higher density of semi-flexible polymers. The total density differences 
between the coexisting phases is however only about 1%. The pairwise inter- 
molecular interactions V/,/(r) (/, J=A,B) are treated as point interactions of 
strength ezijS{r) /<^q. zjj parameterizes the local fluid structure of the under- 
lying microscopic model, as discussed above. The coupling between individual 
chain conformations and the coordination numbers, which results in the spatial 
dependence of the ^-parameter observed in the simulations, is neglected. Fur- 
thermore, we ignore purely entropic contributions (which have been determined 
to be small by Monte Carlo simulations) and do not include orientation depen- 
dent segmental interactions, which will eventually lead to a nematic phase at 
much higher bending energies /. Thus we take the interactions to be 

1^\(PA + <PB - I) ^'Pa 2~ ^ eZAB(PA(PB (9) 

The inverse compressibility C, has been measured in simulations of the athermal 
model; C = 4.1|Q. A Hubbard-Stratonovich transformation rewrites the many 
chain problem in terms of independent chains in external, fluctuating flelds Wa 
and Wb- 

j V[WA,WB.^A,^B]e-^^{-HWA,WB,^A,^B]/kBT) (10) 

where the free energy functional is deflned by 

T[Wa,Wb,^a,^b] 4>A, 2 . , 2 , W^3 c/^ ^ ^ 

[ d^r{WA'i>A + WB^B}-^\liqA[WA]-^lnqB. 
V J I\a j^b 

(j)A = "^^^ =1 — 05 denotes the average A-monomer density and gAiW^] the 
single chain partition function in the external fleld Wa 

QA [Wa] = / ^1 M'Pa [r] exp (^-$o J d^- $aWa^ (12) 
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respectively. The leading contributions to the partition function stem from those 
values 4>A,4'B, WA , of the collective variables which extremize the free energy 
functional, and the mean field approximation amounts to retaining only these 
contributions. The values are determined by: 

-r— =0 Wa = YJ- / d^"^ E{(j)A,(t>B) = C,{4>A + - 1) - eZAA4>A + ^Za^T^ 

d(pA 0(pA J 

n ^ J. (hVSqA 

=0 =^ 9A^^j T (14) 

6wa NaQa owa 

and similar expressions for wb and (f>B- The saddle point integration approx- 
imates the original problem of mutually interacting chains by one of a single 
chain in an external field, which is determined, in turn, by the monomer den- 
sity. Composition fluctuations are ignored, but the coupling between chain 
conformations (e.g. orientations) and the monomer density is retained. The free 
energy of a homogeneous system takes the Flory-Huggins form: 

i'A , iiA , T ^ 1 f/ , r, , \ 72 



In {(j)A)^ T-j In (l - 0A)-:re {{zaa + '2zab + zbb) 't'A~'2 (zab + zbb) (f^A + zbb] 



^okBTV Na ' ' Nb ' ' 2 

(15) 

where we identify the Flory-Huggins parameter x — {zaa + 2zab + zbb)^/'^- 
In the strongly segregated regime, the free energy of a system containing one 
interface is given by: T ^ -t{(j)AZAA + {^-4>A)zBB)^okBTV/2 + akBTL'^. The 
definition of the interfacial tension as the difference of the free energy of a system 
containing an interface and the homogeneous bulk system corresponds literally 
to the measurement of the interfacial tension via the reweighting scheme pO[| in 
the Monte Carlo simulations. As shown in Sec. lie for the symmetric blend, 
these values agree with the measurement of the interfacial tension via the cap- 
illary fluctuation spectrum, so that we can compare the results of the SCF 
calculations and the values extracted from the capillary wave spectrum quanti- 
tatively. 

For the special cases of Gaussian chains jsj, |3p and worm- like polymers [p^ 
^ one can treat the single chain problem in an arbitrary external field in lim- 
iting cases (e.g N — > oo) analytically. For general parameters, however, one 
has to resort to numerical procedures even for these simple models. The BFM 
chains used in the simulations are characterized by structure on different length 
scales. The conformations are rod-like for length smaller than the persistence 
length, which depends on the bending energy /. On intermediate length scales, 
they obey self-avoiding walk statistics, while on the largest scale, the excluded 
volume interactions are screened in the melt, and the chains exhibit Gaussian 
statistics. Since we want to explore dependence on the explicit chain structure, 
we evaluate the single chain partition function via a partial enumeration scheme, 
introduced by Szleifer and coworkers [p^. The method is conceptually straight- 
forward and applicable to arbitrary architecture JlO], |20|] . It can use experimental 
or simulational data as input. Note that no adjustable parameters are involved 
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in the chain structure (such as the statistical segment length in the Gaussian 
model or the bond length and the bending rigidity in the worm-like polymer 
model) and the chain structure is correctly represented on all length scales. Us- 
ing Monte Carlo simulations of the pure melt, we generated 40,960 independent 
polymer conformations for each bending energy. Rotating and translating those 
original conformations, we obtain a sample of 7,864,320 polymer conformations 
for chain length N — 32. (Note only the z coordinates of the chains are em- 
ployed for a flat interface parallel to the xy plane.) For = 64 we use twice as 
many conformations. Within this framework, the A-monomer density (c.f. eq. 

is simply the statistical average of independent A-polymers in the external 
field WA- 





(-Et>A(r«,)) 


J2a=l exP (- WA{ra,i)] 


1 



(16) 



Other single chain quantities are given by corresponding averages over indepen- 
dent chains in the fields wa and wb- 

The set of nonlinear equations is expanded in a Fourier series and solved 
by a Newton-Raphson like method. Convergence is usually reached within 3-6 
steps. The evaluation of the partition function in the external fields poses 
rather high memory demands (several Gbytes). Therefore we employ a CRAY 
T3E, assigning a subset of conformations to each processing element. Typically 
we use 64 or 128 processors in parallel, and the program scales very well with the 
number of processors employed pOt. One needs about 1200 seconds for each set 
of parameters, which is roughly 2 orders of magnitude less than for the detailed 
Monte Carlo simulations. 



3 Comparison between Monte Carlo simula- 
tions and self-consistent field (SCF) calcula- 
tions 

In the following we compare our Monte Carlo simulations to the results of the 
SCF calculations. Both, large length scale thermodynamic properties (e.g. inter- 
facial tension) as well as the local interfacial structure (e.g. orientation of indi- 
vidual bonds) are investigated. The temperature dependence of most quantities 
for symmetric blends (/ = 0) has been studied previousl y ||30| and compared 
to predictions of the Gaussian and worm-like chain model[|l7||7 The results for 
vanishing bending energy compare well to our calculations. 
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3.1 Interfacial tension 

The interfacial tension a between the coexisting phases has an important impact 
on the morphology of the compound material The control of domain size 
and shape is a key to tailoring the application properties of the blend. The 
size of minority droplets often is the smaller, the smaller the interfacial tension 
between the coexisting phases In the strong segregation limit, Helfand 
and Sapsep^ obtained for infinite long Gaussian chains in an incompressible 
blend the analytic expression: 



The interfacial tension a grows upon increasing the bending energy (i.e. the 
statistical segment length) of the semi- flexible component. This behavior is 
presented in Fig. ^ as well as our simulation results and the SCF calculations, 
which take account of the detailed chain architecture. All data exhibit an in- 
crease of the interfacial tension of about 30%. The simulation data and the SCF 
calculations agree nicely on the growth of the interfacial tension upon increasing 
the bending energy. The almost quantitative agreement indicates that our iden- 
tification of the x-parameter yields reasonable results for structural asymmetric 
mixtures. 

However, the Gaussian chain result is about a factor 1.3 higher than the 
simulation data. Recently, Ermoshkin and Semenov^^ calculated corrections 
to the interfacial tension due to effects of finite chain length N . For symmetric 
blends, they found that chain end effects reduce the interfacial tension by a 
factor (1 ~41n2/x-/V) ~ 0.67, which accounts well for the discrepancies between 
the Helfand-Sapse result and the Monte Garlo data. Similar reductions are 
found in numerical SGF calculations [0 |3^ . 

Note that purely entropic, packing contributions to the Flory-Huggins pa- 
rameter Ax are less than 1% of the total x value for / = 1. That is somewhat 
smaller than the uncertainties in identifying the enthalpic contributions of x 
and the accuracy of our interfacial tension measurement in the simulations. 
Therefore, purely entropic effects derived from packing are irrelevant to the in- 
terfacial behavior for the chain lengths, stiffness asymmetries, and temperatures 
investigated in the present study. 

3.2 Monomer density profiles 

Another important characterization of the interface are the density profiles of 
the individual components. Experiments indicate, that entanglements in 
the interfacial zone are of major importance for the mechanical properties of 
the blend. Of course, our chain lengths are too small to observe entanglements, 
however static properties can be extracted from our simulation data. The den- 
sity profiles obtained from the SCF calculation are presented in Fig. H, as well 




(17) 
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as the "apparent" profiles in the Monte Carlo simulation. The width of the ap- 
parent profile in the Monte Carlo simulations is about a factor 1.5 larger than 
the SCF result; this is not unexpected, because capillary fluctuations increase 
the squared width by a term proportional to ln(L)/cr. However, the profiles are 
qualitative similar: both data show a reduction of the total monomer density 
at the center of the interface (the relative reduction is roughly x/2C[0) and 
almost no dependence on the bending energy /. 

The dependence of the interfacial width on the bending energy is shown in 
Fig. The width Wr of the reduced profile is smaller than the apparent 

width Wa, and agrees better with the SCF results. Due to remaining capillary 
wave effects it is an upper bound on the intrinsic width. All profiles presented 
below are obtained by the reduced averaging procedure. The excess energy 
density of the interface can also be used to estimate the intrinsic width. Since 
the relative increase in interfacial area due to capillary fluctuations is of the 
order cr In L/L^, this quantity (as well as the interfacial tension) is not strongly 
affected by fluctuations of the local interfacial position. A tanh-shaped profile 
(j)A = 5(1 + tanh:^) yields in the SCF framework: 



= <Po dz I ^'i^A 2~ ^ <^zab(Pa(PbJ ^ ~ -^We^oX 

(18) 

where we have assumed incompressibility and neglected the finite range of in- 
teractions and any contribution of the intramolecular interactions to the excess 
energy density. Of course, this measure relies crucially on the identification of 
the Flory-Huggins parameter. However, the method is computational very con- 
venient and can be combined ||3^ with the reweighting methods of measuring 
interfacial tensions. It results in values which are between the reduced width 
and the SCF results, which shows again the consistent parameterization of the 
local fluid structure. The Gaussian chain model for iV — > cx) predicts for sym- 
metric blends (/ = 0) a width which is about 20% smaller than the SCF result. 
SCF calculations [|l^ of Gaussian chains with the same long distance behavior 
and which include chain end effects and the finite compressibility agree within 
2% with our results for symmetric, fiexible mixtures. An increase of the chain 
length from TV = 32 to 64 reduces the effective x-parameter by 4% and reduces 
the broadening due to finite chain length effects. The latter effect is stronger, 
such that the width decreases slightly. 

Most notably, the apparent width of the Monte Carlo data, the energetic 
width and the results of the SCF calculations show almost no dependence 
on the bending energy /, whereas the analytic expression obtained by Helfand 
and Sapse 

/ ;,2 ^ 7,2_ 

(19) 

predicts an increase of about 28% due to the variation oi bs- Taking account 
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of the stiffness dependence of the effective Flory-Huggins parameter x, the for- 
mula above predicts an increase of 21%. Quahtatively, calculations for worm-like 
chains [p^ indicate that increasing the bending rigidity results in a reduc- 
tion of the interfacial width compared to the Gaussian chain result. For the 
present combination of parameters, both effects seem to cancel, resulting in an 
interfacial width, which is nearly independent of the bending rigidity. For lower 
incompatibilities, the interfacial width is larger, the Gaussian description on the 
length scale of the interfacial width becomes more appropriate. Therefore the 
difference between the width in the flexible/semi- flexible blend and the width 
in the symmetric flexible mixture increases in accord with the Helfand Sapse 
description. This is confirmed by SCF calculations (c.f. Fig. where we have 
assumed that the effective coordination numbers are temperature independent. 
However, upon increasing the incompatibility further (e > 0.082), one finds 
that an increase of the statistical segment length results in a smaller interfacial 
width of the asymmetric blend in qualitative contrast to the predictions of the 
Gaussian model. 



3.3 Distribution of chain ends and orientations 

The enrichment of chain ends at the center of the interface |^ and at hard 
walls p6|, 48, 49| has attracted abiding interest. Chain ends are important for 



the interdiffusion and healing properties at interfaces between long polymers ||6|. 
They also play an important role for reactions at interfaces. In many experi- 
mental systems, chain ends have slightly different interactions than inner chain 
segments, which might result in a modification of the interfacial properties. On 
the theoretical side, the behavior of chain ends is related to corrections to the 
ground state approximation. Therefore it is a sensitive test for a quantitative 
theoretical description. Chain end effects give rise to large corrections to the 
interfacial width and tension, and they also play an important role for long 
range interactions between interfaces|Q. The distribution of chain ends for 
symmetric blends has been investigated by Monte Carlo simulations!^, and in 
the framework of SCFT for Gaussian chains[^ |l^. In Fig. ^ the simulational 
results and the SCF calculations are presented; both agree almost quantita- 
tively. As in symmetric blends, chain ends are enriched at the center of the 
interface, and this effect goes along with a depletion away from the interface. 
The fact that the depiction zone in the wings shifts outwards with increasing 
chain length, indicates that the length scale of the rearrangement of chain ends 
is the radius of gyration. A-polymers stick their ends into the B-rich phase and 
vice versa. The effect on the semi-flexible chains becomes more pronounced with 
growing stiffness, while the A-polymers are hardly influenced by the stiffness of 
the B-polymers. 

The instantaneous shape of a polymer coil is a prolate ellipsoid [^. Polymers 
orient themselves by putting their ends preferentially at the center of the inter- 
face. This is quantified by the orientational parameter |30| for the end-to-end 
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vector (cf. Fig. 10) 



where the outer index z at the brackets denotes the z coordinate of the midpoint 
of the end-to-end vector R, and the inner indices its Cartesian components. The 
chains ahgn their two long axis parallel to the interface in their majority phases, 
similar to the behavior at a hard wall. The chain orientation of semi-flexible 
polymers increases for growing stiffness, while the flexible A-polymers are not 
affected. The agreement between Monte Carlo simulation and SCF calculations 
is again almost quantitative. In the SCF framework, the orientations of the 
chains in the minority phase is accessible. The polymers align perpendicular 
to the interface, as to reach with one end their corresponding bulk phase. The 
length scale of the ordering increases with the bending energy / and with chain 
length N. 

The orientation of individual bond vectors shows a similar behavior. Bonds 
align parallel to the interface; the effect for the semi- flexible component grows 
with increasing bending energy and its range is largely independent of the chain 
length. The agreement between simulations and SCF calculations is very good. 
The Gaussian chain model cannot predict any nonzero orientation of the bonds. 
The orientation of bonds in our model is, in fact, much smaller than for the 
end-to-end distance Q. 

In contrast to the width of the density profile, the spatial range over which 
the orientation of bonds extends grows upon increasing the bending energy. 
Therefore, the orientational width and the width of the composition profile are 
two independent microscopic length scales. 



4 Conclusions and outlook 

In summary, we have presented extensive simulations of highly incompatible 
polymers with different stiffness. The local structure of the interface has been 
characterized by density profiles of different monomer species and chain ends 
and orientational profiles of whole chains and individual bonds. The interfa- 
cial tension has been measured via analyzing the spectrum of capillary fluc- 
tuations. Using the pair correlation functions of the pure components and a 
random-packing like assumption for the intermolecular contacts between dif- 
ferent species, we have extracted an effective Flory-Huggins parameter, which 
takes account of the stiffness dependence of the structure of the polymeric fluid. 
The effective Flory-Huggins parameter grows upon increasing the stiffness, be- 
cause back folding is less probable and the number of intermolecular contacts 
increases respectively. 

This effective Flory Huggins parameter was then employed in SCF calcula- 
tions, as well as the chain conformations in the pure melt. These calculation 
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incorporate the chain structure on all length scales via a partial enumeration 
scheme; there is no free parameter in describing the chain architecture. Using 
the detailed local structure of the bulk (as obtained by simulations) in the SCF 
calculations, we predict the interfacial properties. 

Monte Carlo results and SCF calculations for the interfacial tension, the 
excess interfacial energy, the redistribution of chain ends and orientations of 
whole chains and individual bonds agree very well provided that the analy- 
sis accounts for capillary fluctuations. However, comparing our results to the 
analytical predictions of the Gaussian chain model for infinite chain length, 
we find qualitative deviations, especially for the dependence of the interfacial 
width on the chain stiffness. This finding might be important for extracting the 
Flory-Huggins parameter from interfacial profiles in highly incompatible poly- 
mer blends. Therefore, our results emphasize that the local structure, both of 
the underlying monomer fluid and of the chain architecture, is important for a 
quantitative description. 

The radius of gyration determines the range of orientation of whole chains 
and the distribution of chain ends. Furthermore, we identify two independent 
microscopic length scales of the interfacial profile; one controls the width of the 
monomer density profile, the other corresponds to the range of orientations. 
This behavior resembles the findings in symmetric blends of worm-like chains 
in the limit kx S> 1 ]l6[ and the behavior of a homopolymer melt at a hard wall 
which is the limiting case for infinite incompatibility. However, in the present 
study this behavior is found for a different model which can be described nei- 
ther by Gaussian nor by worm-like statistics on small length scales. Deviations 
from the Gaussian model occur under rather mild conditions which correspond 
roughly to kx = 0.18 • • • 0.47 in the equivalent worm-like chain model. Further- 
more our self consistent field approach as well as the simulation techniques are 
applicable to arbitrary chain architecture^]. 

Assuming that the chain conformations and the local fluid structure are ap- 
proximately independent of temperature, we have extended the self consistent 
field calculations to other incompatibilities. The results indicate that chain 
architecture becomes important when its length scale is comparable with the 
interfacial width. At very high incompatibility, increasing the stiffness of the 
semi-flexible component results in a decrease of the interfacial width. How- 
ever, the Gaussian chain results and our calculations, which take account of 
the explicit chain architecture on all length scales, agree better for lower in- 
compatibilities, where the interfacial width is much larger than the persistence 
length. 
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Table 1: Single chain conformations and interfacial properties as a function 
of the bending energy /. Interfacial data refer to blends of flexible (/ = 0) 
and semi-flexible (/ as indicated) chains. (6^): mean squared bond length, R^: 
mean squared end-to-end distance, i?^ mean squared radius of gyration, zbb 
effective intermolecular coordination number as measured in simulations of the 
bulk system, {e)/kBT: bulk energy density, eg/kBT: interfacial energy excess 
per unit area, Wa apparent interfacial with for L = 64, Wr reduced interfacial 
width for block size B = 16, We estimated interfacial with from the excess 
energy, Wsd interfacial width in the SCF calculations. 




Figure 1: Single chain structure function S{q) for bending energies / = 0.0,0.5,1.0,1.5, and 2.0 (solid 
lines) . The dashed line represents the Debye function. The arrows, on top of the figure, characterize the 
interfacial length scale (2'u;scf = 6.8). 

Left inset: Mean squared end-to-end distance (circles) and radius of gyration (squares) as a function of 
the bending energy /. Upon increasing /, the chain extension grows about a factor 1.5. Note that the 
ratio Rl/Rg remains close to its Gaussian value 6. 

Right inset: Kratky plot of the single chain structure factor S{q) for / = 0.0,0.5,1.0,1.5, and 2.0 from 

top to bottom. A plateau, characteristic of Gaussian statistics, is only found for / = (flexible chains), 
whereas there are pronounced deviations for semi-flexible chains. 



200 




Figure 2: Scaling plot for the intermolecular pair-correlation function for bending energies / = 0, 1 and 
2. The inset shows the pair-correlation function for the lowest and the highest bending energy studied. 
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Figure 3: Reduced profiles of the coordination numbers of inter- and intrachain contacts for the bending 
energy / = 2 The flexible component is on the left side of the interface. The line presents the apparent 
profile for the AB intermolecular contacts. The dashed line is the random mixing like approximation for 
the AB intermolecular contacts. The inset shows the total number of contacts and number of interchain 
contacts. 




Figure 4: Probability distribution of Fourier components of the local interfacial profile for the 4 smallest 
wave- vectors q and bending energies / = and 2, using the 1/q^ dependence of the variance. The solid 
line represents the expected Gaussian distribution. 

inset: Inverse width of the distribution <7{q) — ^2^2^n2) for all bending energies (circles / = 0.0, squares 

/ — 0.5, triangles / — 1.0, crosses / — 1.5, and stars / — 2.0) and the smallest q values. The 3 smallest 
wave- vectors were used to determine the interfacial tension a. The Helfrich Hamiltonian describes the 
data down to smaller q values for the symmetric blend than for / = 2. The solid line denotes the 
independent estimate for the symmetric blend (/ = 0) from simulations in the semi-grandcanonical 
ensemble. 




Figure 5: Intcrfacial tensions a for all bending energies. Circles represent Monte Carlo estimates obtained 
from the capillary fluctuation spectrum, squares denote the values of the self consistent field calculations, 
and solid lines show the Helfand-Sapse prediction (where we have neglected the architecture dependence 
to the Flory-Huggins parameter i.e. x = 0.265 for all /). 
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Figure 6: A- and B-monomer density, and total density profiles for bending energies / = 0, 1, and 2 
from SCF calculations for chain length N — 32. The apparent profiles in the Monte Carlo simulations 
are shown in the inset). The Monte Carlo profiles are broader than the self consistent field results, but 
in both cases there is almost no dependence on the bending energy /. 




Figure 7: Interfacial width as a function of bending energy / for chain length N ~ 32. Circles denote 
the apparent width Wa in the Monte Carlo simulations {L = 64), squares represent the Monte Carlo 
result Wr from the reduced profiles [B = 16), diamonds show the estimate from the interfacial excess 
energy, filled squares mark the results of the self consistent field calculations, and the solid line shows 
the Helfand-Sapse prediction (x = 0.265 for all /). The inset presents the same data normalized by the 
interfacial width of the symmetric blend. 
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Figure 8: Difference between the width of the / — 0// — 2 blend and the symmetric / = mixture 
in the self consistent field calculations. For e < 0.082 the interfacial width increases upon increasing the 
statistical segment length of the semi-flexible component, whereas at higher incompatibilities the width 
decreases. 
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